Ab initio calculations of the electronic structure of cuprates using large scale cluster 

techniques 



S. Renold, C. Bersier, E. P. StoUQ and P. F. Meier 
Physics Institute, University of Zurich, CH-8057 Zurich, Switzerland 

The local electronic structures of La2Cu04, three members of the Yttrium- family (YBa2Cu306, 
YBa2Cu307, and YBa2Cu408), and to some extent of Nd2Cu04 have been determined using all- 
electron ab-initio cluster calculations for clusters comprising up to thirteen planar copper atoms 
associated with their nearest planar and apical oxygen atoms. Spin-polarized calculations in the 
framework of density functional theory have enabled an estimation of the superexchange couplings 
J. Electric field gradients at the planar copper sites are determined and their dependence on the 
occupation of the various atomic orbitals are investigated in detail. The changes of the electronic 
field gradient and of the occupation of orbitals upon doping are studied and discussed. Furthermore, 
magnetic hyperfine fields are evaluated and disentangled into on-site and transferred contributions, 
and the chemical shifts at the copper nucleus are calculated. In general the results are in good 
agreement with values deduced from experiments except for the value of the chemical shift with 
applied field perpendicular to the Cu02-plane. 
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I. INTRODUCTION 

The discoverjii of high temperature superconductivity 
has initiated great experimental and theoretical efforts 
which aimed at a detailed understanding of these ma- 
terials. Nevertheless, a generally accepted theory which 
explains at least the most important properties of the 
high temperature superconductors could not yet be pre- 
sented which is possibly related to the complex electronic 
structures of these materials. 

Therefore, very early on ab-initio methods have been 
employed to determine the electronic properties of these 
materials. Mostly, band-structure techniques have been 
employed which are reviewed in Ref. 0]. In our con- 
trasting approach to electronic structure calculations, we 
use cluster techniques in the framework of spin-polarized 
density functional theory with localized basis functions 
which are especially well suited for the calculation of lo- 
cal properties. 

In particular, the cluster method has been successfully 
applied to the determination of charge and spin-density 
distributions in the cuprate plane and to a reasonably ac- 
curate evaluation of electric field gradients and magnetic 
hyperfine fields for planar copper and oxygen nuclei in 
La2Cu04 and YBa2Cu307 using clusters comprising five 
planar copper and their nearest neighboring planar and 
apical oxygen atoms^i^. 

In this work, these calculations are extended to in- 
clude clusters representative of La2Cu04, YBa2Cu306, 
YBa2Cu307, YBa2Cu408, and to some extent of 
Nd2Cu04. The development of both hard- and software 
has led to a significant improvement of the quality of 
our calculations by use of larger clusters comprising up 
to thirteen planar copper atoms. This allows a careful 
study of the convergence of the calculated local proper- 
ties with respect to the cluster size. 

The paper is organized as follows: We will first outline 



the cluster technique and introduce the used clusters in 
Sec.ini Sec.|lll]is devoted to a discussion of the spin dis- 
tribution and contains estimations of exchange couplings 
in the cuprates. In Sec. lIVI we report on the calculation of 
electric field gradients (EFG) at the planar copper sites 
in the different cuprates. The convergence of the results 
with respect to the cluster size is demonstrated and the 
various contributions of the molecular orbitals (MO) and 
atomic orbitals (AO) to the EFG are elucidated. The 
calculated distribution of charges and holes reveals that 
the 3o?322_^2 AO of the copper is not fully occupied. The 
doping dependence of these distributions and of the EFG 
are evaluated and compared to nuclear quadrupolc res- 
onance data. Sec. |V] contains a careful examination of 
hyperfine parameters and their disentanglement into on- 
site and transferred contributions. Calculations of chem- 
ical shieldings are presented in Sec. IVII and in particular, 
the role of the reference substance for magnetic shifts of 
the copper is discussed. The paper is terminated with 
the summary and conclusions in Sec. IVIII 



II. THE CLUSTER TECHNIQUE 

The idea of the cluster technique is to select a con- 
tiguous region out of the solid and to treat the elec- 
trons therein with standard many-body theories. This 
so-called core region is surrounded by a shell of basis- 
free pseudopotentials which prevent the electrons to be 
attracted by positive point charges and provide smooth 
boundary conditions. The core and the boundary re- 
gions are embedded in a large lattice of background point 
charges to ensure a very good approximation for the 
Madelung potential. In the case of the two undoped par- 
ent compounds, La2Cu04 and YBa2Cu306, the values 
chosen for the background point charges are based on 
the formal valence of the constituents. In clusters rep- 
resentative of YBa2Cu307 and YBa2Cu408 the formal 
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valences had to be slightly modified to reach charge neu- 
trality. 

The calculations were performed in the framework of 
spin-polarized density functional theory which provides a 
good trade-off between accuracy and computational cost. 
(Hartree-Fock theory, which completely neglects correla- 
tion effects, yields poor results. In contrast, configura- 
tion interaction methods, which correctly include both 
exchange and correlation effects in their Hamiltonian, 
require enormous computer resources and are currently 
prohibitive for the cluster sizes considered here.) For 
the representation of the exchange and correlation func- 
tionals, the potentials of Becked and Lee, Yang, and 
Parr—, have been used. All atoms in the core region of 
the clusters employ the standard triple zeta basis sets (6- 
311G). For the determination of the ground state of the 
many-electron system and the evaluation of the observed 
quantities, the GaussianOS quantum chemistry software 
package^ was used. 

It is desirable to have as many atoms as possible in the 
core region, but the available computer resources and the 
convergence of the self-consistent field procedure are the 
limiting factors in this respect. As already mentioned, 
these limits have been pushed further since our first use 
of the cluster technique (see Ref. [1]) which now allows to 
use larger clusters including up to nearly 1000 electrons. 

In this work we present results for clusters compris- 
ing 5, 9, and 13 planar copper atoms together with 
their nearest planar and apical oxygen atoms for clusters 
representative of La2Cu04, YBa2Cu306, YBa2Cu307, 
YBa2Cu408, and Nd2Cu04. In Table [Til the constitutive 
properties of all the used clusters for every substance are 
listed. The layout of the Cu02 plane for the three clus- 
ter sizes is displayed in Fig. 2. The lattice constants and 
the positions of the atoms in the unit cells for the three 
clusters of the Y-family were chosen according to ex- 
perimental structure determinations and are taken from 
Refs. [ilillllllil- The buckling of the planar oxygen 
atoms is not shown in Fig. 2. For La2Cu04, the calcu- 
lations were performed for the tetragonal structure with 
lattice constants a = b = 3.77 A, c = 13.18 A and atomic 
positions according to Ref. (loj . 



III. SPIN DISTRIBUTION 

Within the spin-polarized formalism, the spin multi- 
plicity M of a cluster is free parameter of the calculation. 
In a simple ionic picture, the planar copper and oxygen 
atoms have a valence of +2 and —2, respectively This 
leads to a 3dP configuration for the copper atom with a 
total spin of one half whereas the oxygen valence gives 
zero total spin. This suggests two choices of spin mul- 
tiplicities with a physical significance, a "ferromagnetic" 
spin multiplicity with all spins parallel and an "antifer- 
romagnctic" spin multiplicity with two neighboring spins 
being antiparallcl. However, other spin multiplicities are 
also possible leading to spin alignments which can be 



substance 


name 


N 


P E 


B 




CU5O26 


31 


42 395 


533 


La2Cu04 


Cug042 


51 


62 663 


897 




CU13062 


75 


78 991 


1313 




CU5021 


26 


194 345 


468 


YBaaCusOe 


CugOss 


42 


62 473 


780 




CU13049 


62 


86 841 


1144 




CU5021 


26 


37 345 


468 


YBa2Cu307 


CugOss 


42 


62 473 


780 




CU13049 


62 


86 841 


1144 




CU5021 


26 


37 345 


468 


YBa2Cu408 


CU9033 


42 


53 473 


780 




CU13049 


62 


73 841 


1144 



TABLE I: Compilation of the used clusters with their con- 
stitutive properties. N: number of atoms with a full basis 
set, P: number of atoms with (basis-free) pseudopotentials, 
E: number of electrons in the core region, B: number of basis 
functions. 
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FIG. 1: (color online). The layout of the CUO2 planes in all 
of the used clusters, (a) with 5 planar copper atoms, (b) with 
9 planar copper atoms, and (c) with 13 planar copper atoms. 
The apical oxygens are not shown. The empty triangles de- 
note Cu^^ ions simulated by a pseudopotential. 



viewed as superpositions of an antiferromagnetic and a 
ferromagnetic spin state. We anticipate that the calcu- 
lated total energy is always lowest for M that corresponds 
to an antiferromagnetic alignment. In Table HITl the cho- 
sen spin multiplicities M are tabulated for each cluster 
size. The antiferromagnetic (ferromagnetic) spin multi- 
plicities are in bold (normal) face. The superpositions 
are in parentheses. 

cluster multiplicities 



CU5O26 (2), 4, 6 

CU9O42 2, (4), (6), (8), 10 

CU13O62 (4), 6, (8), (10), (12), 14 

TABLE II: Chosen spin multiplicities for the different clus- 
ters. Numbers in bold face denote multiplicities leading to 
an "antiferromagnetic" (4, 2, and 6) and a "ferromagnetic" 
(6, 10, and 14) spin alignment. Note that the possible spin 
multiplicities are determined solely by the number of planar 
copper atoms and are therefore independent on the specific 
material. 
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In Fig.[2]thc spin density along the planar Cu-0 bonds 
is drawn for the large Cui3 cluster representative of 
YBa2Cu306 in the case of M = 14 corresponding to 
the ferromagnetic spin alignment (upper panel) and in 
the case of Af = 6 corresponding to the antiferromag- 
netic spin alignment (lower panel). The total energy for 
the latter spin multiplicity is 2.9 eV lower than that of 
the former. The double humps at the copper positions 
originate from the approximately singly occupied 3dx2-y2 
atomic orbital. The spins of neighboring coppers are in- 
deed parallel for the ferromagnetic case and antiparallel 
for the antiferromagnetic case thus confirming the sim- 
ple physical picture given above. It is important to note 
that these spin density distributions are obtained also in 
the smaller clusters and in the other substances consid- 
ered. A close inspection of the values at the copper sites 
(upper panel) shows a difference between the three inner 
coppers and the two coppers at the borders. This is due 
to the fact that the latter have only one nearest-neighbor 
(NN) Cu ion which implies a single transferred hyperfine 
field {B > 0) whereas the formers have 4 NN leading to 
AB. This will be further discussed in detail in Sec. IV Al 
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FIG. 2: Plot of the spin density along the planar Cu-0 bonds 
in the Cuis cluster for YBa2Cu306. The upper (lower) panel 
shows the ferromagnetic (antiferromagnetic) spin configura- 
tion with M = 14 (M = 6). 

Furthermore, the ground state energy of the antiferro- 
magnetic spin alignment is consistently lower than that 
of the ferromagnetic spin alignment. This important fea- 
ture of our cluster method can be exploited to investigate 
on the exchange couplings in the different materials. 

The Heisenberg Hamiltonian is given by 



(1) 



where the summation is restricted to pairs of indices 
with the corresponding sites being nearest neighbored 
(and with j > i to avoid double counting). 

For a given material and cluster size, we can obtain 
ground state energies Em and ground state wavefunc- 
tions V'M(r, (t) = (^Af(r)XA^ ("■) for any of the allowed 



multiplicities M (see Tablc|TTT|. This is, however, not suf- 
ficient for a determination of J since ( xm I ' I Xm ) 
is in general not easy to determine. However, for a two 
spin state, the Heisenberg model has an exact solution 
with (X3| Si • S2 I X3 ) = 3/4 and ( Xi I Si • S2 | Xi ) = -1/4. 
This leads to the following idea: 

We replace the expectation values of the spin opera- 
tors in the Heisenberg Hamiltonian by the product of the 
MuUiken spin densities at the respective lattice sites. 



:XM|S,-S,|xM)=pf (CuOpf (Cu,), 



(2) 



with the requirement that in the triplet state of the two- 
spin system 



ap3(Cu,)p^(Cu,) = 1 



(3) 



with a reduction factor a. 

(Note that the name "spin density" might be mislead- 
ing. It is in fact not a spin density, but an integrated 
spin density, i.e. a spin. Nevertheless, we persevere with 
this term which is commonly used in quantum chem- 
istry.) With clusters comprising only two planar copper 
atoms, we determined /9s(Cui) = 0.67 for both copper 
sites which yields a = 2.23. This reduction of the spin 
density from its value in an ionic model (p^(Cui) < 1) is 
due to the fact that also the planar oxygen atoms carry 
a small amount of spin density. It has nothing to do 
with the reduction of the mean magnetic moment in the 
Heisenberg model by quantum fluctuations. We thus ar- 
rive at the following modified Heisenberg-typc equation 
for the determination of J: 



E 



M 



-a J 



{ i-3 ) 



M 



(Cu,)pf (Cu, 



-aJS. (4) 



In Fig. [3] the ground state energy is plotted versus S 
for the different spin multiplicities M in the case of the 
large Cuia cluster representative of La2Cu04. A straight 
line with slope aJ = 350 meV can be fitted. 
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FIG. 3: Ground state energy versus E (see Eq. |4j for all 
spin multiplicities M of the Cuia cluster representative of 
YBa2Cu306. 

The results of this fit are given in Fig. 4. It is observed 
that the exchange couplings are fairly independent on the 
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size of the cluster. Furthermore, they are all in the same 
range of about 150 meV. It is only the optimally doped 
YBa2Cu307 that has a slightly lower exchange coupling 
which is mainly due to the buckling of the planar oxy- 
gens. The two antifcrromagnetic substances La2Cu04 
and YBa2Cu3 06 and the underdoped YBa2Cu408 have 
similar values of J. 

The calculated antifcrromagnetic couplings arc in good 
agreement with datapi^^^^ and also with theoretical val- 
uesi^ obtained from smaller clusters but with more so- 
phisticated ab-initio methods than used here. 

We therefore conclude that the antifcrromagnetic ex- 
change interactions between two copper neighbors are an 
intrinsic property of the Cu02 plane and depend only 
weakly on the specific material of experiments. 



analysis of the contributions to the copper EFG in terms 
of molecular orbitals is presented. We then postulate 
in subsection IIVDI an approximate relationship between 
the copper EFG values and the partial atomic MuUiken 
charge populations of the orbitals 3(i^2_j^2 and Sd^^^^r^. 
Although approximate within about 10 % only, this pro- 
vides insight into the EFG differences observed in dif- 
ferent compounds as well as their changes upon electron 
or hole doping. The latter relationship is discussed in 
subsection IIVFI 

In subsection llV GI thc theoretical results are compared 
to experimental data and an extended discussion about 
the interpretation of the results is presented. 

B. Dependence of the calculated Cu EFG on the 
cluster size 
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FIG. 4: Exchange couplings fitted as described in the text. 



In Fig. 5 we present the main component Vzz of the 
EFG obtained for the copper in the center of clusters 
comprising = 5, 9, and 13 planar copper atoms sim- 
ulating the La2Cu04, YBa2Cu408, YBa2Cu306, and 
Nd2Cu04 compounds. These values have been calcu- 
lated for spin multiplicity M = 4, 2, and 6, respectively, 
which correspond to an antifcrromagnetic alignment of 
the spins and yield the lowest total energy. The values 
for N = 9 and 13 are nearly equal while those for N = 
5 are about 10 % to 20 % smaller due to finite size ef- 
fects. We note that the EFG values for the YBa2Cu307 
compounds exhibit the same size dependence with val- 
ues for YBa2Cu307 between the those for La2Cu04 and 
YBa2Cu306 (see subsection IIV Fp . 



IV. ELECTRIC FIELD GRADIENTS 
A. Outline 

Electric field gradients (EFG) are extremely sensitive 
to the (non-spherical) charge distributions around the nu- 
cleus of interest. In the past twenty years a large quan- 
tity of nuclear quadrupole resonance data has been accu- 
mulated for high-temperature superconducting cuprate 
compounds. Of particular interest are the EFGs at the 
planar copper site and their changes upon doping since 
they reflect the local charge distribution and provide in- 
sight into the population of the different atomic orbitals. 

We have already published calculated values for the 
EFGs in La2Cu04, YBa2Cu307, and Nd2Cu04 obtained 
for clusters with five copper atoms {N = 5) in the plane 
(see e.g. Refs. [IQ)- In this paper we add results for 
larger clusters {N = 9 and 13) and report in addition 
results for YBa2Cu306 and YBa2Cu408. 

In subsection IIV Bl we first investigate on the influence 
of finite size effects. Next (subsection IIV C|) a detailed 
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FIG. 5: EFG value V,, 
atoms in the clusters. 



as a function of the number of Cu 



In addition, we have calculated the EFG value for a 
La2Cu04 cluster with an even number {N = 12) of cop- 
per atoms without spin polarization (singlet state with 
AI = 1). As seen in Fig. 5 this value is very close to 
those calculated for M = 9 and 13. 

We further note that Vzz is given in atomic units (1 a.u. 
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corresponds to 9.71525 x 10^"'^ V m ^). A comparison to 
experiments will be made in subsection IIV G[ 



C. Detailed analysis of contributions 

In this subsection we first analyze the various contri- 
butions to the copper EFG in detail. It will be shown 
that the values of the EFG are not solely determined by 
the occupancies of the id,^2_y2 atomic orbital as is often 
assumed in simplistic models. 

The charge distribution is determined by the occu- 
pied MOs. The m*^ MO is represented as a linear 
combinatio»ii 

n n nx 

</.„,(f) = -^"i^^- rk) = E E^™'^^>fc(^^- 

K=l K=l k=l 

(5) 

ofriK atomic basis functions i?A',fe centered at the nuclear 
sites K = 1, . . . , n, and the c^'*^ are the MO coefficients. 
We assume in the following that the target nucleus Kt 
is at Rkt ~ 0- The contribution of the MO (pm to the 
EFG at Kt is given by the matrix clement 

n n riK 
K=l L=l fe=l 1=1 



X ( BK.k{r^ Rk)\ I BL.i{r~ Rl) )■ (6) 

This matrix element contains contributions from basis 
functions centered at two nuclear sites K and L. Thus 
we can identify three types of contributions: (i) on-site 
terms from basis functions centered at the target nucleus 
{K = L = Kt, contribution I), (ii) mixed on-site off- 
site contributions (II), and (iii) purely off-site terms with 
K ^ Kt and L ^ Kt (HI). 

In addition, there is a contribution coming from all 
nuclear charges Zk and all point charges mentioned in 
Sec. II which wc denote by 

where Wy refers to nuclei in the core region of the cluster 
and Wf^'^ to the point charges outside this core. The 
charges of the bare nuclei at sites K ^ Kt in the core arc 
screened by the matrix elements of contribution III with 
K = L ^ Kt- Therefore the combined contributions 
from III and the nuclei (Wy ) are small. The partitioning 
of the contributions to the EFG tensor Vij thus reads 

V^, = 'V,, + "V,, + '"v., + W,, + lyP^ (8) 



More details about these regional partitions can be found 
in Ref. [H. 

In Table IIIII the contributions to the EFG component 
Vzz for the copper in the center of the Cuia clusters rep- 
resentative for Nd2Cu04, YBa2Cu307; and La2Cu04 are 
collected. As expected, the values from region III and W 
as well as those from point charges outside the core region 
are small. However, mixed on-site off-site contributions 
(region II) are substantial. They are mainly transferred 
via the on-site 3(i3j,2_^2 orbital which has a non- negligible 
overlap with the 2pz orbitals of the four surrounding pla- 
nar oxygens. The on-site terms (region I) contain values 
which come from terms in the evaluation of matrix ele- 
ments ^ where one basis function is s like and the other 
d like. Their contribution is denoted by R in Table V. 
Also the contributions from the Cu p type orbitals are 
substantial due to the large ( ) values. The occupan- 
cies of the three d orbitals with t2g symmetry are close 
to 2 so that their combined contribution to Vzz is small. 
As expected, the distinguished AO is the M^^-y^ with a 
polarization of 70 % accompanied with smaller but still 
significant polarization in the orbital Zd^^'^-r^- 

D. Model 

As pointed out above, a rigorous explanation for the 
variation of the copper EFG values in different cuprates is 
very complicated. A simple explanation, however, which 
concentrates on the two most relevant orbitals, is given 
in this subsection. 

Previously, we have investigated in Refs. [T9ll20j the 
changes of the copper EFG values that occur upon dop- 
ing in La2Cu04 and in Nd2Cu04. We found that it is 
sufficient to concentrate on those MOs which contain par- 
tially occupied 3c?^2_j,2 and 8^322 _^2 AOs at the target 
nucleus. The EFG is then given by 

Vzz^ 8+^^ [iV3;,2_^2(r"^)3^2_^2 - N^2_y2{r~'^)^2_y2] 

(9) 

where S is the contribution from all other orbitals and 
regional partitions. The values for (r~^) calculated for 
the different Sd atomic orbitals are almost identical (e.g. 
(r"3)3j,2_r2 = 8.004 a.u. and {r-^)^2_y2 = 8.042 a.u.) 
and therefore are replaced by the average < >= 
8.020 a.u.. We further note that the partial occupation 
numbers N of the AOs are very similar to the partial 
MuUiken populations pc which are gathered in Table IIVI 
In particular, the differences 

Ad = Pc(3c?3z2-r2) " Pc(3da;2_y2) (10) 

and 

■^3z2-r2 - Nx2_y2 (11) 

are very similar. In Fig.[S]we therefore plot the calculated 
copper EFG component Vzz versus A^. The straight line 
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TABLE III: Contributions to the EFG component for 
the copper in the center of the Cuia clusters representative 
for Nd2Cu04, YBa2Cu307 and La2Cu04. 





Nd2Cu04 


YBa2Cu307 


La2Cu04 




-0.033 


0.021 


0.019 


I 


R 


0.755 


0.640 


0.474 


P 


-1.431 


-1.157 


-1.250 


d^2_y2 


-6.378 


-6.344 


-6.263 




-0.072 


-0.141 


-0.090 


d322_^2 


8.309 


8.630 


8.736 


II 


s 


-0.003 


-0.004 


0.004 


p 


0.013 


0.014 


0.022 


d^2_y2 


-0.021 


-0.020 


-0.019 


d3;i/~!~d3^2~t~djy2 


0.001 


0.001 


0.002 




-0.585 


-0.515 


-3.382 


III + Nuclei 


-0.015 


-0.061 


-0.053 


Total 


0.538 


1.066 


1.202 



TABLE IV: Partial occupation numbers N and partial Mul- 
liken populations pc for the frontier orbitals 3d^2„j,2 and 

3(^322 _^2. 

Nd2Cu04 YBaaCugOT La2Cu04 



"FT 



1.379 
1.430 
1.820 
1.865 



1.375 
1.423 
1.898 
1.914 



1.363 
1.406 
1.910 
1.922 



The corresponding values for the four units surrounding 
the central one are 

Pi(Cu) = 1.131, pi(Op) - -1.629,pi(0,) = -1.944. 

(14) 

while those of units at the boundary of the core region 
differ at most by one percent due to finite size effects. 
We define 



corresponds to 

V,, = -2.74 + 7.63Arf (12) 
which emphasizes the quality of the model. 
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0.42 0.44 0.46 , 0.48 0.50 0.52 



Po(3) = p(Cii) + 2p(Op) + 2p(0,)+6 (15) 

where 6 accounts for the charges of the two La ions in 
the imit cell. We obtain po(3) = —0.012 that is close 
to zero which emphasizes the suitability of these clusters 
to represent the local conditions of atoms in the crystal. 
This point cannot be overstressed since it also means 
that there is now also confidence in the partial MuUiken 
populations of the individual orbitals which have been 
studied in detail elsewhere [2^ . 

The partial MuUiken populations of those AO which 
differ from 2 by more than 0.02 are given in Table IVl and 
are visualized in Fig. 7 where the area in yellow (gray on 
top) accounts for an intrinsic hole of 1.5 missing electrons 
which is compensated by the occupancy of the 4s orbital. 
All attempts to deduce the distribution of holes which do 
not take into account the 4s orbital are very questionable 
and misleading. 

4 I \ 



FIG. 6: (color online). Theoretical EFG component 
Vzz versus Ad=pc(3d3^2_r2)— Pc(3d^2_y2) for La2Cu04 
(black square), YBa2Cu307 (blue triangle up), YBa2Cu408 
(blue triangle down), YBa2Cu306 (blue triangle right), 
Sr2Cu02Cl2 (from Ref [IJ]) (green circle), and Nd2Cu04 (red 
diamond). The straight line is a fit V-^z = 7.63 • Ad — 2.74. 
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E. Charge and hole distribution and qualitative 
discussion of bonding 

The atomic MuUiken charges p for the central copper 
and neighboring oxygen atom in La2Cu04 as calculated 
for the Cui3 with spin multiplicity M = 6 are 

Po(Cu) - 1.143, po(Op) = -1.635, po(Oa) = -1-945. 

(13) 



FIG. 7: (color online). Partial MuUiken populations of the 
atomic orbitals (tinted in darker gray (blue)) for La2Cu04: 
Black symbols denote 4s (circles), ?>d^n_yn (diamonds), and 
3d322_^2 (triangles up), 2 x p{0^(2pcr)) for the planar oxy- 
gens (squares), and 2 x p(0a(2p2)) for the apical oxygens 
(triangles down). The light gray (yellow) tinted areas denote 
the missing charge against the simple ionic model. 
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TABLE V: Partial MuUiken populations pc of the AOs in the 
central unit of La2Cu04. 



AOs 




Pc 


1.406 1.922 0.491 1.662 1.946 



The ionic picture which assigns charges of +3(La), 
+2(Sr), +2(Cu), and —2(0), respectively is a reason- 
able approximation for the out-of-plane La- and apical 
Oa-atoms. The copper and the planar oxygen atoms, 
however, are rather covalently bound and the relation 
between charge and hole transfer is not trivial. A Mul- 
liken charge population analysis [23| attributes a charge 
of 1.16 to the copper and —1.64 to the planar oxygens 
but 40% of the hole is transferred to the oxygens leaving 
60% on the copper. Since a hole transfer is accompanied 
by an electron charge transfer in the opposite direction 
it is concluded that of the total of 0.84 electrons trans- 
ferred from the oxygen to the copper, 0.40 electrons are 
accounted for the hole transfer. The remaining 0.44 elec- 
trons almost correspond exactly to the MuUiken 4s or- 
bital population. The discrepancy of 0.05 electrons can 
be attributed to secondary interactions which also involve 
the 3(^322 _^2 and the 2pz apical oxygen orbitals. The sub- 
stantial occupation of the Cu 4s orbital is responsible for 
the existence of the hyperfine field transferred from the 
four next nearest copper atoms which is revealed by NMR 
experiments. 

F. Doping dependence 

The measured copper quadrupole frequencies gen- 
erally increase with hole doping. In particular, in 
La2-2:Sr2,Cu04, the slope of this increase is reported 
i n^3^24,25 |.jjg range of 0.56a; to 0.7a;. It is of ut- 

most interest to deduce the redistribution of charges that 
occurs upon doping from these data. This is, however, 
a complex task because the doped holes will also change 
the lattice parameters and atomic positions in the unit 
cell which in turn influence the EFG. 

Ab-initio calculations of doping-induced charge re- 
distribution have been reported by Ambrosch-Draxl et 
al. for the compound HgBa2Cu04+i-. They em- 
ployed the full-potential linearized augmented plane- 
wave method and used a series of supercclls containing 
one excess oxygen atom. In principle, total-energy and 
atomic-force calculations can also be performed for small 
clusters. They are not feasible, however, for large clus- 
ters. We therefore report in the following on investiga- 
tions of doping dependence with lattice parameters and 
atomic positions kept fixed at values corresponding to the 
undoped compound La2Cu04. The results are therefore 
only reliable for small doping, i.e. they should be consid- 
ered to describe the (linear) slope of changes around the 
undoped material. We will turn back to changes in the 
atomic positions at the end of this subsection. 

The method of cluster calculations would also enable 



the study of changes in the local electronic structure that 
occur by replacing e.g. a trivalent La by a bivalent Sr. 
These inhomogeneous changes, however, are not the sub- 
ject of the present investigation where we simulate the 
long range effects of doping (delocalized holes or elec- 
trons) by two different approaches. First we introduce 
additional point charges at the periphery of the cluster 
to place an electric field across the cluster to move charges 
toward or away from the atoms of interest in the cluster 
center. In this so called "peripheral charges method", 
the added system of charges has no physical interpreta- 
tion except that it can be continuously altered so that 
charge can be progressively directed to or extracted from 
important ions of interest in the cluster. The cluster and 
the system of charges in total keep the same number of 
electrons and spins but, using for instance a MuUiken 
population analysis approach, the charge of the cluster 
center can be progressively changed in a manner expected 
by doping. In the second approach of simulating doping 
we have added or removed two electrons from the cluster 
and repeated the calculation. Changing the number of 
electrons by an even number allows us to keep the spin 
state. 
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d(3)=p(3)-pJ3) 



FIG. 8: V^;, versus "doping level" d(3) = p(3) - po(3) for 
La2Cu04. 

To discuss the EFGs for doped La2Cu04 we define 

p(3) = p(Cu) + 2p(0p) + 2p(0a) + 2p(La), (16) 

where the MuUiken charges refer to the central unit and 
p(La) — 3. 

For the undoped Cui3 cluster with spin multiplicity 
M=Q we obtain po(3) = —0.012 as discussed above. In 
Fig. [H] we plotted the EFG component Vzz for the central 
copper versus the doping level d(3) = p{2>) — po(3)- Note 
that negative values of d would apply to electron-doped 
materials which for the present case of La2Cu04 is of 
course of no experimental relevance. The two circles at 
d — —0.09 and d = 0.11 have been obtained by adding 
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1.44 



1.40 




1.94 



1.90 



d(3)=p(3)-p^(3) 



FIG. 9: Pc{Sd^2 _y2) (solid curve) and Pc{Sd3^2_j.2) (dashed 
curve) versus "doping level" d{3) — p(3) — po(3) for La2Cu04. 



and subtracting two electrons, respectively. The trian- 
gles denote results obtained with the peripheral charge 
method [l^. All values have been calculated for multi- 
phcity M=6. 

The slope of the increase of Vzz at d — is 1.12. 
An analysis of the populations of the frontier orbitals is 
shown in Fig.O For small d{3) but with increasing (hole) 
doping Pc(3d3z2_r2) is increased but reaches a maximum 
at d{3) = —0.056 while Pc{3dx2-y2) is decreased and has 
a minimum at 0.066. Using Eqs. dH), ^ and the 
dashed line in Fig. [8] shows the calculated Vzz- 

An analysis of the charges in the occupancies upon 
doping exhibits that an additional extrinsic hole goes to 
15 % to the 3d^2_y2 AO, to 2 x 18.7% = 37.4% to the 
two Op 2pa- AO, to 7 % to the Sd^z^-r^ AO and to 2 x 
18.7% = 37.4% to the two Oa 2pz AO, while the 4s 
remains practically constant. 

It is instructive to compare these results with those 
obtained in Ref. [1^ for IIgBa2Cu04+5. Although the 
calculational procedures arc quite different the essential 
conclusions arc the same. We first define 



was reported whereas our value for La2Cu304 is 15 %. 
This is compensated by a smaller decrease of 20 % for 
2 X pc{2p„) compared to 37 %. These differences are, 
however, of not too much relevance since the assignment 
of the charge to AO in covalent bonds is anyhow some- 
what arbitrary. The turn-over of Vzz in Fig. 8 and in the 
occupancies in Fig. 9 may well be an artefact of the not 
appropriately adjusted changes in the atomic positions. 
It should be pointed out, however, that in the calcula- 
tions of Ambrosch-Draxl et al. [1^ where these positions 
have been adjusted, the initially linear changes stop at 
a doping concentration d{i) = 0.22 and remain constant 
at higher (i(3) values which indicates a saturation of the 
planar hole content at 0.12. 

The essential question now is whether these theoreti- 
cally predicted redistributions of charges can be corrobo- 
rated by experimental facts. The slope of increase of Vzz 
of 1.1 is too large compared to the data. We may explain 
this partially with arguing that the lattice parameter a 
shrinks upon doping. While it is not feasible to deter- 
mine the ground-state energy of all atoms with cluster 
calculations, it is straightforward to study the changes if 
a single parameter is varied and the results are reliable 
since relative changes are involved. We have previously 
reported^! on the changes that occur when the lattice pa- 
rameter a for a cluster representative of La2Cu04 with 
5 atoms is varied and found that when a shrinks by 1 % 
the EFG value Vzz is reduced by ~ 10.7 %. Assuming 
a reduction of a by 4 % upon doping, according to data 
presented in Ref. [2I], we get a reduced slope of Vzz of 
1.12 — 0.43 = 0.69, which is (very) close to the experi- 
mental value. Notice, however, that this argumentation 
neglects the variation of Vzz due to small changes of the 
positions of the atoms in the c-direction, in particular 
those of the apex oxygen and of the La ions. 



G. Comparison with NQR experiments 

For a nuclear spin 3/2, the connection between NQR 
frequencies ^^vq and the main component Vzz of field 
gradients in the case of axial symmetry is given hy^: 



p{2) = p{U,2_y2) + 2p{2p„) (17) 

as the charge of the planar orbitals and correspondingly 
d{2) = p(2) — po{2) as the derivation from the undoped 
case. We get 

p{2) ^ 0.52 p{3) (18) 

which is to be compared with p{2) = 0.55 p(3) in Ref. 
[2^. Thus in both cases the removal of an electron by a 
dopant atom in the intra-layer induces only half a hole 
in the Cu02 plane. The calculated changes of the occu- 
pancies of the individual orbitals are somewhat different. 
For HgBa2Cu04+5 a decrease of 35 % for pc (3^3,2 _,y2) 



63 Ai^(±3/2 ^ ±1/2) e ^^QVzz 

with ^^Q being the nuclear quadrupole moment. Unfor- 
tunately, directly measured values for ^^Q are not avail- 
able. The commonly used value of^^Q = — 0.211 ± 
0.004 b was obtained some time ago by Sternheimer— 
by interpreting excitation spectra with the Hartree-Fock 
approximation. 

In Fig. (TU] we plot our theoretical values for the 
quadrupolar frequency ^^vq obtained from the calculated 
values Vzz with the above mentioned value for ^^Q versus 
the calculated difference of the MuUiken populations A^. 
The straight line corresponds to 
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0.42 0.44 0.46 0.48 0.50 0.52 



FIG. 10: Experimental (full symbols) and calculated (open 
symbols) quadrupole frequencies ^^uq as a function of the cal- 
culated difference of the MuUiken populations A^. The sym- 
bols denote: La2Cu04 (square), YBa2Cu307 (triangle up), 
YBa2Cu408 (triangle down), YBa2Cu306 (triangle right), 
Sr2Cu02Cl2 (from Ref [2l[) (circle), and Nd2Cu04 (dia- 
mond) . 



^^UQ =a + PAd, a = -73 MHz, /3 = 200 MHz. (20) 

In the same figure wc include data obtained from 
the various materials plotting them at the correspond- 
ing theoretical values A^. The experimental values for 
^'^lyg are very accurate but there is of course an uncer- 
tainty in the calculated A^. Furthermore, we expect the 
theoretical values for Vzz to be less reliable for small 
Vzz due to the cancellations of the various contributions 
shown in IIV CI The more reliable calculations for the 
compounds with large Vzz (La2Cu04, YBa2Cu307 and 
YBa2Cu408) systematically predict quadrupole frequen- 
cies that are about 15 % lower than the experimental 
ones. The dashed straight line corresponds to 

= ci + /3Arf-/(a + /3Arf). (21) 

with / = 1.2 = 1/0.83. This disagreement between cal- 
culations and experiments could be due to a systematic 
error in the theoretical determination of Vzz or due to a 
higher value of ^^Q than assumed or due to both. 

It should also be noted that the NQR frequencies ^^i^q 
depend on temperature as has been reported in detail 
by Matsumura et al. ^l| for La2Cu04. According to 
the cluster calculations with variable lattice parameter a 
(see IV F) the general increase of ^^Vzz with temperature 
in the paramagnetic region is mainly due to the lattice 
expansion. 

Owing to the orthorhombic structure, the EFG at 
the planar Cu(2) is not axially symmetric in the com- 
pounds YBa2Cu307 and YBa2Cu408 but shows a slight 



anisotropy since Vxx 7^ ^yy The calculated anisotropy 
parameters t] = \Vxx — Vyy\/\Vzz\ are r] = 0.047 and 0.035, 
respectively. 

The experimentally observed anisotropies have been 
reviewed by Brinkmann^. For the planar copper nu- 
clei, they are somewhat smaller than our theoretical val- 
ues. Of more importance, however, are the large rj values 
(slightly below 1) measured for the chain copper, Cu(l), 
which are quite unexpected since Cu(l) is not at a crys- 
tallographic position that would imply 77 = 1 by sym- 
metry arguments. An evaluation of the EFG at Cu(l) 
by cluster methods requires that at least the two near- 
est neighboring planar copper atoms in the Cu02-planes 
above and below the Cu(l) site are considered. We have 
previously performed^i^ such calculations for a CU3O12 
cluster representative of YBa2Cu307. The calculated^ 
EFG values for Cu(l) are Vxx = 0.601, Vyy = -0.603 
and Vzz = —0.002 a.u. which should be compared with 
the experimental ones Vxx = 0.767, Vyy = —0.773 and 
Vzz = 0.006 a.u. as obtained from the measured fre- 
quencies with the above mentioned quadrupole moment 
^^Q. The theoretical values produce an asymmetry pa- 
rameter in complete agreement with the data whereas 
the absolute value are again about 20 % smaller than the 
measured ones. 

Recently, Kanigel and Keren^ have reported 
NMR measurements for a series of fully enriched 
(Cao.iLao.9)(Bai.65Lao.35)Cu30y powder samples where 
doping can vary across the full range from the very 
underdoped to the extreme overdoped. They determined 
the nuclear quadrupole frequency from the four peaks 
of the powder spectra and obtained a convex curve of 
i^Q versus doping level which looks very similar to the 
one (triangles) shown in Fig. [H They used then also 
a simulation program to account for asymmetric peaks 
which provides an increase of vq with doping in the 
underdoped region but a saturation at the overdoped 
site. 



V. MAGNETIC HYPERFINE FIELDS 

A. Theoretical determination 

The hyperfine spin Hamiltonians for copper and oxy- 
gen nuclei in the Cu02 planes of the cuprates are given 
by: 

s^i/f-^ = I, • A • S, + h-M-Sj (22) 

and 

''^f^= E I'-C-S.- (23) 

jeNN 

For copper, the hyperfine interaction contains an 
anisotropic on-site term. A, and a transferred term, B, 
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from the four nearest copper neighbors, which is usually 
taken to be isotropic, since it is assumed to consist only 
of a contact term. We will, however, retain the tensorial 
notation of the transferred field, since we will also be able 
to calculate a dipolar part of the transferred interaction. 
For oxygen, the hyperfine interaction C is with the two 
nearest copper neighbors. 

We have already demonstrated in Refs. [113] that for 
La2Cu04 and YBa2Cu307 it is justified to neglect trans- 
ferred interactions from further distant copper neighbors, 
both in the case of copper and oxygen. 

On a first-principles level, the hyperfine interactions 
are basically well known, they consist of an isotropic hy- 
perfine density, D, a dipolar interaction, T*^ , and a spin- 
orbit interaction term. The core polarization is given by 
the spin density at the nuclear site R and is evaluated as 
follows: 

D(R) = ^(pT(R)-pi(R)). (24) 

If it mainly originates from singly occupied s-clcctrons, 
it is called Fermi contact. In contrast, if it is due to a 
doubly occupied s-statc which is polarized through the 
spin of other electrons at the same atom or on remote 
atoms, it is called core polarization. Since in our case, 
both contributions are present, we prefer to call D just 
the isotropic hyperfine density. 

The dipolar interaction is evaluated as 



T,,(R) 



dMp^r) - pHr))A^(r) (25) 



with 



A^(r) = 



3{n-R,){rj^Rj)-S,j\r-Il\' 



r- 



R|5 



(26) 



The estimation of the spin-orbit interaction will be dis- 
cussed later. The above equations yield the total hy- 
perfine fields and still have to be split into on-site and 
transferred contributions. To exemplify this splitting we 
first focus on the clusters with maximal multiplicity. In 
our clusters for YBa2Cu306 we find copper atoms with 
no nearest neighbors (in a small Cui cluster), with one 
(the corner copper in the largest Cuis cluster), with two 
(e.g. the corner copper in the intermediate Cug cluster), 
with three (the edge copper in the Cuia cluster) and four 
nearest neighbors. We plot in Fig. [TT] the value of the 
isotropic hyperfine density against the number N of near- 
est copper neighbors and find a linear dependence of D 
on the number of nearest copper neighbors. The value of 
D when there is no nearest copper neighbor present, is 
then the on-site contribution, and the slope of the fitted 
straight line in Fig.[TT]is the transferred contribution per 
copper neighbor according to the ansatz 



The numerical values for these two contributions are 
ttiso = —2.00 a,g^ and biso = 0.52 a^'^. In the very same 
way, one can obtain values for the on-site and transferred 
part of the dipolar interaction. The z-components are 

II —3 II —3 

'^dip — —3.38 a^ and b'^-^^ = 0.06 a^ . It should be em- 
phasized that in our quantum-chemical calculations the 
on-site and the transferred terms arc highly connected 
and the linear dependence of D extends from zero NN 
up to four NN. 




FIG. 11: Plot of the isotropic hyperfine density recorded 
at different copper sites in clusters of various sizes for 
YBa2Cu306 against the number of nearest copper neighbors 
at that specific site. 

By a slight generalization of the ansatz ([?7|) it is also 
possible to include results from clusters with lower mul- 
tiplicities in the determination of on-site and hyperfine 
fields. Wc write 



D{CUi) aisoPs (Cu,) + ft.o /'^(Cuj) (28) 

and similarly 

r-(Cu,0 = a^,pP.(Cu,) + /'^(C^^)- (29) 

For the hyperfine fields at the oxygen site, we make a 
completely analogous ansatz: 



and 



^(O,) =7^.o J2 '^-(C^j) 

j£NN 



(30) 



(31) 



D 



(27) 



The connection between the fitting parameters a 
through 7 and the actual hyperfine parameters is then 
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given by scaling with the expected spin density in the in- 
finite cluster for which we find a good estimate in the cen- 
ter of the largest Cuia clusters with ferromagnetic spin 
multiplicity for each substance, i.e. Oiso = Ps{Gu)aiso 
etc.. 

The quality of the ansatz can be estimated by suit- 
ably chosen plots as e.g. a plot of D{Cui)/ ps{Cui) 
against X]jgArArPs(Cuj)/ps(CuO. (see Fig. [T2|). From 
the straight line, the fitting parameters aiso and /^i^o are 
determined and scaled with the expected spin density in 
the infinite cluster - as noted above ~ to get the hyper- 
fine parameters Uiso = —1.94 a^'^ and biso = 0.77 a^"^. 
The corresponding values for the dipolar hyperfine cou- 
pling are a||.p = -3.55 a^^ and bl- = 0.08 a^^. These 
values are very similar to the ones obtained with the sim- 
pler ansatz which shows that the two different kinds of 
ansatz yield effectively the same results. Because in the 
second ansatz more clusters are included we quote in the 
following these results. 



2 I 1 1 1 1 1 r 
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FIG. 12: Plot of D{Cui)/ps{Cui) against 

Jl^jg^jv Ps(Cuj)/ps(Cui) for YBa2Cu306. Closed symbols 
originate from the cluster with ferromagnetic multiplicity 
M = 14 and open symbols denote results for antiferromag- 
netic multiplicity M = 6. Circles are for values where Cui 
is in the center of the cluster whereas points for off-center 
coppers are plotted with triangles. 

In Table IVII all hyperfine parameters as determined 
using the ansatz II can be found. We note that for the 
oxygen c^jp refers to the direction along the bond between 
the two NN coppers, c^^^ is perpendicular to the bond but 
still in the Cu02 plane, while c^f^ denotes the direction 
perpendicular to the plane. 

The contribution to the hyperfine fields that originate 
from spin-orbit coupling a"^ are expected to be small in 
the case of oxygen. For copper, however, they are of 
the same order of magnitude as aiso and a^ip. We are 
not in a position to determine them at the same level of 
accuracy as Oj^o and adip and are thus forced to rely on 
a reasonable estimate. In the frame of an atomic picture 
with a single missing electron in the 3dx2-y2 orbital the 



La2Cu04 YBaaCusOe YBaaCuaOy YBa2Cu408 





-1.94 


-2.09 


-2.03 


-2.08 




0.77 


0.57 


0.50 


0.55 


a". 

dip 


-3.55 


-3.38 


-3.40 


-3.38 


^dip 


0.08 


0.05 


0.06 


0.06 


C-iso 


0.64 


0.63 


0.62 


0.60 


^dip 


0.40 


0.40 


0.39 


0.42 


^dip 


-0.20 


-0.21 


-0.21 


-0.22 


^dip 


-0.20 


-0.19 


-0.18 


-0.20 



TABLE VI; Theoretical values for the hyperfine parameters 
at the planar copper and oxygen sites in the four substances 
determined with the more sophisticated ansatz II using Cuia 
clusters with ferromagnetic and antiferromagnetic spin ar- 
rangement. All values are given in units of a^''. 



dipolar and spin orbit hyperfine interactions are given by 
(Ref. [H) 

aLp = -y^"^). alo = -yfc^"^> anda^-^ = -yfc(r-3). 

(A value of fc = —0.044 was estimated in Ref. [36[). In 

a molecule or solid where the missing electron spends 

some time on the oxygen ligands these expressions have 

to be modified. A simple modification is to replace the 
II .... 

expression for a^^^^ by multiplying it with 2 — N^2_y2 . 

Using the values {r'^) = 8.042 (Sec. IVD) and N^2_y2 = 

1.363 (Table V) we obtain for La2Cu04 ajj^p = -2.93 
which is reasonably close to the directly determined value 
of —3.55 since actually the spin density should be used 
instead of the charge density. 

For the estimation of the spin-orbit contributions we 
thus assume that the relations ((5^ still hold in the cluster 

such that alo = 15.5A:a|-p and = 2.75ka\^p. The 
resulting values are given in Table IVIII 

In Table IVIIII we collect the calculated values for the 
total hyperfine parameters expressed in terms of densities 
and also in terms of interaction energies which we denote 
by capital letters. The latter are defined by ^^Af^^ = 
hjefi^^^a'tot ^^rd similarly for the ^^B and ^^C and (as 
the notation) indicates depend on the particular isotope. 



La2Cu04 YBaaCugOe 


YBaaCugOr 


YBa2Cu408 


alo 2.43 2.32 


2.27 


2.28 


aj-o 0.43 0.41 


0.40 


0.40 



TABLE VIL Estimation for spin-orbit contributions to the 
hyperfine fields at the copper site in the various substances. 
All values are given in atomic units. 



B. Comparison with experiments 

Although experiments cannot determine on-site and 
transferred hyperfine fields separately it is possible to 
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U.Ul 
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U.OO 


0.D2 


O.OD 


O.Oi 




0.73 


0.52 


0.47 


0.52 


Hot 


1.04 


1.03 


1.01 


1.02 


<hot 


0.44 


0.42 


0.41 


0.38 


Ctot 


0.44 


0.44 


0.44 


0.40 


63 /ill 

A-tot 
63 4± 

63 oil 

^tot 
63 D± 

t>tot 
17^11 

17,^± 

17,^c 


-1.79 
0.15 
u.ou 
0.43 
0.31 
0.13 
0.13 


-1.84 
0.01 

U.OD 

0.30 
0.31 
0.13 
0.13 


-1.85 
0.04 
n 

U.OO 
0.27 
0.30 
0.12 
0.13 


-1.86 
0.01 

u.ou 

0.30 
0.30 
0.11 
0.12 



TABLE VIII: Contributions to the hyperfine fields at the cop- 
per site (a and 6) and at the oxygen site (c) in the various 
compounds in &^ {a,b and c) and in ^eV ^^B and 

"C). ^'^A/a and '^^B/b = 0.5844 /ieVa| and ^'^C/c = 0.2988 
/ieVa^ . 



extract various combinations of on-site and transferred 
fields using different experimental set-ups. 

One constraint is the relation a^^^ + ^bf^^ = which ex- 
plains that the NMR spin shift measured with the field 
in c-direction does not change below the superconduc- 
tion transition temperature Tc in contrast to the spin 
shift measured with the field perpendicular to the c-axis. 
We postpone a comparison with our theoretical values to 

sec. ivrc] 

A second experimental determination of a combina- 
tion of on-site and transferred hyperfine fields is made 
possible through measurements of the NMR resonance 
frequency, ^^vli of the copper nuclei in the pure par- 
ent compounds La2Cu04 and YBa2Cu406 which are in 
the antiferromagnetic state. This determines the local 
magnetic field ^^jHioc = iir^^VL which is commonly ex- 
pressed as the corresponding hyperfine field in units of 

an effective electronic magnetic moment ~ 

as Hioc = \ai^^ - '^b^J^if. 

For La2Cu04 the theoretical value for A = ja^^ — 4&^j| 
is 2.66 whereas the measured^ frequency of 93.85 MHz 
corresponds to A = 2.01. For YBa2Cu306 we get 
A = 1.90 in good agreement with the experimenlj^^ 
(89.89 MHz) which leads to A = 1.92. 

Sometimes the measured anisotropy between the cop- 
per spin-lattice relaxation times ^'^Ti n/^^Ti j^ has also 
been used to extract information about the hyperfine cou- 
pling constants. This, however, is only possible for the 
two extreme conditions of totally antiferromagnetic cor- 
relations (which is never achieved for the doped samples) 
or of no correlations, which would require measurements 
at very high temperatures. 

In Table IIXI we compare our values for the hyperfine 
interaction energies with those published by several au- 
thors. An inspection shows that the deviations are not 





a) 


h) 


c) 


d) 


e) 


63 4 II 

A-tot 

63 

A-tot 
63 n 


-1.85 
0.04 
0.29 


-1.86 
0.23 
0.34 


-1.76 
-0.10 
0.41 


-1.61 
0.29 
0.40 


-0.94 
0.17 
0.23 



TABLE IX: Contributions to the hyperfine fields at the copper 
site in YBa2Cu307. a) this work, 6) Ref. [s^, c) Ref. (3^, d) 
Ref. HOl, and e) Ref. gl[. 

large but sufRcicntly strong to render further interpre- 
tations questionable. In particular, until more reliable 
values for the influence of the spin-orbit coupling on the 
hyperfine fields are known, it is prohibitive to make more 
precise statements. 

VI. CHEMICAL SHIELDINGS AND 
PARAMAGNETIC FIELD MODIFICATIONS 

A. General remarks 

Very early after the discovery of the high tempera- 
ture superconductors, numerous measurements of Knight 
shifts at various nuclei have been performed. For opti- 
mally doped YBa2Cu307, the Knight shift ^^K^ of the 
planar copper for the field applied perpendicular to the 
c-axis is temperature independent above Tc and drops 
below Tc with decreasing temperature to ^^K^{T = 0). 
The behavior above Tc is to be expected for a tempera- 
ture independent Pauli spin susceptibility. The reduction 
of the Knight shift below Tc was explained by the forma- 
tion of Cooper pairs which are - due to their vanishing 
total spin - not available for polarization by a magnetic 
field and therefore do not contribute to the Knight shift. 
At zero temperatures, all charge carriers were assumed 
to be bound in Cooper pairs and the remaining Knight 
shift ^^K^{T ~ 0) was attributed to the temperature and 
doping independent chemical shift. The temperature de- 
pendence of ^^K^{T) in the superconducting state de- 
pends on the symmetry of the pairing state and most 
NMR experiments were better explained by d-wave pair- 
ing. (For underdoped materials, the decline of the Knight 
shift ^^K^{T) with lowering temperatures sets in already 
at temperatures above Tc which has attributed to the 
opening of a spin pseudogap.) 

The copper Knight shift ^^XH (T) with the apphed field 
along the crystallographic c axis, however, is constant 
over the whole temperature range, i.e. it is completely 
unaffected by the superconducting transition. Therefore, 
it was concluded that all of the measured Knight shifts 
are of chemical origin. The vanishing spin part of the 
Knight shifts was then explained by an accidental can- 
cellation of the on-site and transferred hyperfine fields. 
As already mentioned in Sec. |V] our calculations of on- 
site and transferred hyperfine fields in substances of the 
La and Y families do not support this cancellation. In 
this section we will, in addition, give further evidence 
that the above sketched explanation of the Knight shifts 
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in cuprates needs careful revision based on first-principles 
calculations of chemical shifts at the planar copper nuclei 
in La2Cu04, YBa2Cu306, and YBa2Cu307. Before we 
report on the results of these calculations, we point out 
that copper Knight shift measurements have long been 
misinterpreted due to a wrong assumption on the mag- 
netic properties of the reference substance. As a remedy 
we introduce a new term, the paramagnetic field modifi- 
cation. 



B. The role of the reference substance 

Knight shift measurements are in fact measurements of 
differences in resonance frequencies of a nuclear species 
fc, in a target substance (t) and in a reference substance 
(r). 



^K(t-r) 



V(r) 



(33) 



For copper Knight shift measurements the most often 
used reference substance is the monovalent CuCl. The 
Knight shifts are made up of two parts, a temperature 
independent chemical shift, '^i^^,, and a spin shift, ^K^. In 
this section we consider the contribution of the chemical 
shift. 

For a theoretical determination of chemical shifts we 
need to know the chemical shieldings, <t, in the target and 
the reference substance. The connection to the measured 
chemical shift, ^K'f^{t — r), is given by 



^Kt{t-r) 



1 + ^a"{r) 



'=CT**(r)- (34) 



or, with the separation of '^cr" into diamagnetic (d) and 
paramagnetic {p) parts of the shieldings, by 

^Kt{t -r)= ^a^^{r) + 'a;{r) - ^a^^it) - 'a;{t). (35) 

As we have shown in Ref. [i^ the differences between 
^^CT^'(r) and ^^cr^J {t) are negligible and we can therefore 
write 



^^K2it - r) = ^V;*(r) - «V;*(i) 



(36) 



It has long been assumed, that the paramagnetic contri- 
bution of the shielding in the target substance, CuCL is 
small. We have shown, however, (i) by a direct quan- 
tum chemical calculation and (ii) by referring to mea- 
surements of the chemical shieldings using atomic beam 
techniques, that ^"^(7^ (CuCl) = 1500 ppm. In view of 
typical chemical shieldings at copper nuclei in cuprates 
(see Sec. IVI C|) this contribution is sizeable and cannot 
be neglected. 

The quantities of interest in Eq. ([36|) are, of course, 
not the chemical shifts, ^^K"{t — r), but rather the con- 
tributions of the target, i.e. — ^•^(Tp(t). To avoid misinter- 
pretations we find it most convenient to introduce here a 
new quantity, the paramagnetic field modification: 
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-'^a;{t)^'^K^^-r)-'^^a;(r). (37) 



It is important to note that this paramagnetic field mod- 
ification ^^K'l{t) is independent of the reference sub- 
stance. 



C. Theoretical determination of chemical shifts 
and comparison to experiment 

For a detailed description about the determination of 
chemical shifts in the framework of the cluster technique, 
we refer the reader to Refs. [i^lZsj . Here we just mention 
the important steps. 

Calculations have been performed for La2Cu04, 
YBa2Cu406, and YBa2Cu307 using clusters with five 
copper atoms in the plane. Larger clusters have not been 
employed systematically since the determination of the 
shielding constants with the large basis sets that are re- 
quired for accurate calculations are extremely time con- 
suming. Test calculations in Cug clusters for La2Cu04 
have shown, however, that the results do not change upon 
enlarging the clusters. 

In Fig. [T3] theoretical results for ®'^A'|[ (^^i^^) are dis- 
played with solid bars in the left (right) panel. The 
dotted bars denote results reported from various exper- 
iments. It is observed that for the applied field in the 
plane (right panel of Fig. [T^ the theoretical values are 
in general slightly lower than the values obtained from 
experiments. For fields parallel to the c axis, theory pre- 
dicts values for ^'^K^j^ (left panel) that are only about half 
as large as the experimental results. We find, however, 
that the paramagnetic field modifications at the planar 
copper sites hardly depend on the specific cuprate. 

We have, at present, no explanation for the discrep- 
ancy between theory and experiment. We would like to 
point out, however, that it was shown in Refs. (solfsil that 
the measured Knight shifts below Tc are field dependent 
and drop when reducing the applied magnetic field H. 
The theoretical calculations, of course, are in the limit 
of _ff — > 0. Furthermore, it is also possible that impu- 
rities induce a finite density of states at T = 0, as was 
proposed in Ref. [31 . Both of the above two ideas imply 
that the measured Knight shifts at T = are not entirely 
of chemical origin but also have contributions from spin 
degrees of freedom. 

The temperature independence of the copper Knight 
shift when measured with the field in c-direction is most 
easily explained by an incidentend cancellation of the on- 
site and transferred hyperfine fields, i.e. af^j -I- ^b^^^ ~ 
0. In Table X we present the values calculated for 
1 4- al^f/Abfgf for the four substances under considera- 
tion. The values are close to zero but differ among the 
various compounds. For La2Cu04, a small positive value 
of 0.1 is obtained, for YBa2Cu307, however, we obtain 
—0.4. This difference is from the theoretical point of view 
easily explained by the fact that the transferred hyper- 
fine field in YBa2Cu307 is smaller than in La2Cu04 
due to the buckling of the planar oxygen atoms. It is 
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FIG. 13: Paramagnetic field modifications at the planar cop- 
per site in the considered substances both for fields along the 

c axis, ^^i^i, (left panel) and for fields in the plane, ^^K^^ 
(right panel). Results obtained theoretically (experimentally) 
are displayed with solid (dashed) bars. The experimental data 
are taken from Ref. Q (La2-a:Sr3;Cu04), from Refs. [islliel 
(YBa2Cu06) and from Refs. (YBa2Cu07-«). 



1 - 




La2Cu04 


0.105 


YBa2Cu306 


-0.278 


YBaaCugOT 


-0.402 


YBazCuiOg 


-0.292 



TABLE X: Theoretical relations for the total hyperfine fields 
at the central copper atom. 



evident that the calculations for one compound may fail 
to give the cancellations necessary for the easy explana- 
tion of the temperature independence of ^^A'". To ex- 
plain the behavior of ^•^/•sTU both in La2Cu04 and in the 
Y-compounds already requires a double coincidence. In 
addition, measurements on the electron doped material 
PrLaCeCu04 by Zheng et al^ also exhibit a tempera- 
ture independent ^^/C'l. In view of the differences in the 
lattice parameters in all three substances it is extremely 
intriguing that these cancellations of on-site and trans- 
ferred hyperfine fields which are basically determined by 
chemistry, occur. 



VII. SUMMARY AND CONCLUSIONS 

We have performed large-scale ab-initio cluster compu- 
tations of cuprates in order to determine the local elec- 
tronic structure. The convergence of these local prop- 
erties with respect to the cluster size is very good. An 
analysis of the charge and spin distribution in terms of 
contributions from the various MO and AO reveals distin- 
guished features in all compounds under consideration. 
First, the copper 3dx2-y2 AO is occupied by about 1.4 
and the •id^z2_^2 AO by about 1.9 electrons. The oxygen 
2pcr AO contains roughly 1.65 electrons. This implies a 
total of 1.4 intrinsic holes per unit which is compensated 
by 0.4 — 0.5 electrons in the copper 4s AO. These partial 
occupancies of the non-spherical AO mainly determine 
the EFG values. The 4s AO is involved in the transferred 
hyperfine field. Good agreement between the calculated 
and measured copper EFG is found. The EFG values es- 
sentially depend on the differences between the MuUiken 
populations of the Sd^^^-r^ f^nd the 3d^2_y2 AO. 

Simulating doping by two different methods shows that 
all these occupancies smoothly change and the general 
trends of changing copper EFG with doping level are re- 
produced. The removal of one electron by a dopant atom 
in the intra-layer induces only half a hole in the Cu02 
plane. The other half is in out-of-plane orbitals. 

Spin-polarized calculations with various spin multiplic- 
ities enabled the determination of the antiferromagnetic 
exchange coupling and the various hyperfine fields. The 
contribution of the spin-orbit coupling, however, has only 
been approximately determined. The calculated total hy- 
perfine fields are in rough agreement with those deduced 



from experiments. The values for the sum a^^^ + 4b[gf 
small but differ considerably among the substances un- 
der consideration. This in sharp contrast to the require- 
ments set by the temperature independence of the copper 
Knight-shift ^^K^^ observed in very different cuprates. 

We conclude that the out-of-plane orbital 8^3-2 and 
the 4s orbital play a more important role than commonly 
assumed. 



are 
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